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SUMMARY 


Preliminary estimates are prepared of the weighting functions of 
the infrared bands of methane, ammonia, and hydrogen for the Jovian 
atmosphere. Together with the infrared spectral emission observations 
of Jupiter, these weighting functions are used to obtain a rough estimate 
of the vertical temperature profile in the Jovian atmosphere. 



1. INTRODUCTION 

The primary objectives of the current research program are: 

(1) To study the meteorology of the planet Jupiter, with particular 
emphasis on the deduction of the atmospheric temperature profile through 
analysis of available infrared and microwave observations of the planet's 
emission spectrum. 

(2) To develop and apply inversion techniques for inferring atmo- 
spheric temperature and constituent profiles of planetary atmospheres. 

During the first quarter, work was started on the first objective 
listed above. The available observations of the Jovian emission spec- 
trum having the required spectral resolution consist of the infrared 
observations of Gillett et al. (1969) and the microwave observations of 
Law and Staelin (1968) and Wrixon et al. (1971). To invert these obser- 
vations and deduce the temperature structure requires a knowledge of 
the absorption (transmission) characteristics of the Jovian atmospheric 
gases that influence the observed spectrum. These are methane and hydro- 
gen — which absorb in the infrared — and ammonia — which absorbs in both 
the infrared and microwave portions of the spectrum. Much of the first 
quarter was spent in developing the transmission models for these gases. 
These transmission models have been used in preliminary calculations to 
obtain a picture of what the weighting functions look like. (The con- 
cept of weighting function is explained in the next section; briefly. 
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the weighting function indicates where in the atmosphere the observed 
radiation originates.) These weighting functions, together with the 
observations were used to obtain a preliminary and rough estimation of 
the Jovian atmospheric temperature profile. 

This work is reviewed in detail in the following section. The last 
section of this report summarizes research plans for the next quarterly 
period. 



2. DISCUSSION 


2.1 Weighting Function Concept 

A quick insight into the inversion process can be obtained from an 
understanding of the weighting function concept. Thus, in this section 
we shall first discuss the weighting function concept in general, then 
discuss specific transmission models for the Jovian gases of interest, 
and then apply these transmission models to approximate Jovian atmo- 
spheric conditions to obtain the Jovian weighting functions. From the 
Jovian weighting functions and the emission spectrum observations, we 
shall obtain a rough estimate of the temperature profile. 

The emitted radiation intensity of a planetary atmosphere at wave- 
number v can be written as 

00 

I(v) = B[v, 0(z g ) ] T(v, z g ) + J B [ v, 0(z) ] [ £ (V?Z;> ] dz (1) 

z 

s 

where B is the Planck function, Q is temperature, z is altitude, T(v,z) 
is the transmittance from the altitude z to the top of the atmosphere, 
and the subscript s refers to the planet's surface. In the presence of 
a cloud cover that acts as a black body, the cloud top altitude replaces 
the planet's surface in the equation. For strongly absorbing gases, 

T(v, z ) may be zero, and only the second term on the R.H.S. of (1) — 
the atmospheric contribution — need be considered. If there is no 
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absorption at all by the atmosphere, then the second term on the R.H.S. 

of (1) equals 0, T(v, z g ) = 1, and the emitted intensity is equal to 

the Planck function at 0(z ), from which the planet's surface tempera- 

s 

ture can easily be retrieved. Such regions of the spectrum with little 
or no absorption are called window regions. 

From the second term on the R.H.S. of (1) it can be seen that the 
emitted intensity as observed at the top of the atmosphere consists of 
contributions from levels within the atmosphere. These contributions 

consist of the Planck function at a level weighted by the vertical 

dT 
dz 

called the weighting function. 

Physical reasoning can be used to demonstrate that the weighting 
function for a strongly absorbing gas will have a peak somewhere between 
the surface and the top of the atmosphere. At the top of the atmosphere, 
the contribution to the outgoing intensity is negligible because the 
density of the emitting (absorbing) gas is low. At the base of the 
atmosphere, the contribution to the outgoing intensity is small since 
the overlying atmosphere absorbs most of the radiation emitted by the 
near surface layer. Thus, the weighting function must have a maximum 
somewhere between the surface and top of the atmosphere. 

The level at which the weighting function peaks is the center of 
the atmospheric layer from which most of the emitted intensity originates 
We show below that, for a gray absorbing atmosphere with exponentially 


as a function of z is 


dT 


gradient of the transmittance, ^ . Thus, 


I 
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decreasing density of absorbing gas with altitude, this peak occurs at 
the level where the transmittance is equal to e \ 


For a gray atmosphere the transmittance from the level z to the top 
of the atmosphere is given by 


T = exp j - a I p dz 

S 


00 

/ 


( 2 ) 


where a is the absorption coefficient, and p is the density of the absorb- 

§ 

ing gas. If the absorbing gas has a constant mixing ratio 


w = p /p 
g 


(3) 


where p is the atmospheric density. Equation (2) can be transformed with 
the aid of the hydrostatic equation 


dp 

dz j 


(4) 


to 


/? (z) 

T = exp J -k / dp ] = exp 




(5) 


where k = 


QW 


The weighting function is 

I 

*1 = -ke' kp -2 = ncrk e' kp 

dz ke dz Pg ke 


( 6 ) 
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To determine at what level the weighting function peaks, we differen- 
tiate (6) with respect to z, and obtain 



The weighting function is a maximum where 


g k) 

d_ 

dz 




0 or where 


1 

P 


dfi = 
dz 


P g k 


For an isothermal atmosphere, 


( 8 ) 


1 d£ = 1 d£ = _ 1 
p dz p dz H 

where H is the scale height. Thus, (8) is satisfied when 


( 9 ) 


H = P 8 k (10) 

or when k = 1/p. When k = 1/p, the transmittance T = e The opacity 
or optical depth, x, is defined as (-in T) and is equal to 
unity at the level of the peak weighting function. 

We shall later use the level of peak weighting function to obtain 
an estimate of the height from which most of the emitted radiation orig- 
inates, and from the magnitude of the emitted radiation obtain an esti- 
mate of the temperature at this height. Since the level of peak weighting 
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function varies with the absorption coefficient, observations at wave- 
lengths with different absorption coefficients will yield information 
on the temperature at different heights. Such a usage of the weighting 
function represents a rudimentary form of inversion. More sophisticated 
inversion techniques are discussed by Conrath (1968), Chahine (1968), 

Wark and Fleming (1966), and Smith et al. (1970), among others. In future 
work, we plan to utilize an iterative inversion technique of the type 
described by Chahine (1968). Thus, the work reported here represents 
the simplest approach to the problem^that still yields useful results. 

It also has the virtue of yielding insight into the inversion problem. 

The mathematical derivation above applies to a gray atmosphere. To 
determine weighting functions for a non-gray atmosphere, and to invert 
the emission spectrum of a non-gray atmosphere requires a knowledge of 
the real transmission characteristics of the absorbing gases. We discuss 
below transmission models and associated weighting functions for CH^, 

NH^, and on Jupiter. 

2.2 Transmittance Models and Weighting Functions 
2.2.1 7. 7q CH/ j Band 

For the infrared absorption band of CH^ at 7.7n we follow 
McClatchey et al. (1970). McClatchey assumes that the transmission of a 
20 cm ^ interval can be represented as a function of an absorption 
coefficient k^, a path length interval, AL q , and the pressure p, as 
follows: 

V v) " 1 [\ % => a ] < u > 
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From laboratory transmittance data and from line by line transmittance 
calculations degraded to 20 cm ^ resolution, McClatchey derives empir- 
ically the function f. For CH^, he finds n = 0.75. 


The empirical transmittance function is represented in McClatchey 
et al. (1970) by a graph relating transmittance to wavelength for a path 
length of one kilometer of air at standard temperature and pressure (STP). 
The transmittance scale also has a scaling factor which permits reading 
of transmittances for other path lengths. From McClatchey 's . empirical 
transmittance function, we have derived an analytical representation of 
the CH^ transmittance as follows. 


A vertical increment of pressure corrected path length of CH^ can 
be written as 


n. 


CH. 


dL = 


n 


.75 


n STP \ P STP 


dz 


( 12 ) 


where n is the atmospheric number density, and n CH is the CH^ number 
density. 

Letting c = (n^ H ^/n) , and assuming a constant CH^ mixing ratio with 
altitude, we obtain for the pressure corrected CH^ path length of the 
atmosphere above height z 


L = 


.75 

‘STP r STP z 


l n 


.75 a 

p dz 


(13) 
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With the help of the hydrostatic equation in the form 

dp = -n [_i g dz 


( 14 ) 


where |i is the mean molecular mass of the atmosphere, we obtain 


L = 


.75 


/ 


.75 . 

P dp 


n STP P STP ° 


Letting A = 


n STP P STP 


.75 


jig, we obtain, upon integration, 
: 1.75 


L = 


1.75A 


(15) 


(16) 


We shall assume that the transmission follows the law 

T v = exp [-(k v L) m ] (17) 

where is an absorption coefficient for the interval Av centered at v, 
and m is to be determined from the empirical transmittance relationship 
of McClatchey et al. This form is similar to the strong line approxi- 
mation of the Goody random model, except that in the strong line approxi- 
mation the pressure dependence is p\ not p’^, and m = 0.5. 


By plotting in(-in T) versus in X, where X is kilometers of air at 
STP, from McClatchey's empirical transmittance scale, and measuring the 
slope of the resulting plot we .obtain the value of m = 0.56. Values of 
k are determined from the graph presented in McClatchey et al. These 
values are shown in Table 1. 
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TABLE 1 

ABSORPTION COEFFICIENTS FOR METHANE FOR USE IN EMPIRICAL 
TRANSMISSION FORMULA T v = exp f -(l^L) ,56 1 


v(cm _1 ) X(p) k v^ cm STP ^ 


1140 

8.77 

3.1 x 

10 -4 

1160 

8.62 

5.4 x 

10" 3 

1180 

8.47 

5.4 x 

00 

1 

o 

i-H 

1200 

8.33 

6.2 x 

io' 3 

1220 

8.20 

1.7 x 

io‘ 2 

1240 

8.06 

1. 1 x 

io- 1 

1260 

7.94 

6.2 x 

io -1 

1280 

7.81 

9.4 x 

t- 1 

o 

1 

V- 4 

1300 

7.69 

1.4 


1320 

7.58 

2.8 x 

io" 1 

1340 

7.46 

3.1 x 

io" 1 

1360 

7.35 

2.4 x 

h- 4 

o 

1 

1380 

7.25 

1.7 x 

io" 3 
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To obtain the weighting function, we differentiate (17) with respect 
to height, making use of (12) and the ideal gas law. 


u ,. tN -.44 1.75 

dT _ - 56 < k v L) C P 

n _ Potd'^ k 0 
STP *STP 


exp 


[-<V>' 5 ‘] 


(18) 


where k is Boltzmann's constant, and 6 is temperature. 


Using (18) we have computed the weighting functions for the Jovian 
atmosphere. The following assumptions are made: 


_4 

c = n / n = 8.3 x 10 (after Hogan et al., 1969) 

lh 4 

6 = 168°K 
m = 3.26 
0/m = 51.5 

Here 0 is the mean temperature of the atmosphere above the Jovian 
cloud layer. Together, with the mean molecular weight, it determines 
the scale height, 16.5 km. The cloud-top pressure is assumed to be either 
2 atm or 4 atm and the height scale is referenced to the cloud-top. 

The computed weighting functions are shown in Figure 1 . The large 
vertical range in the peaks of the weighting functions suggests that 
temperature information can be obtained for a layer of the atmosphere 
about 100 km thick above the clouds. However, overlap of absorption 
due to other bands (for example, NH^ at X > 8.3|i) may cut down this 
range in practice. Also, noteworthy are the relatively large widths 
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Figure 1. Weighting functions in the Jovian atmosphere for 

selected wavenumbers of the 7.7ji CH 4 absorption band. 
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of the weighting functions, of the order of 40 km at — equal to one- 
half of peak value. This is due to the relatively large scale height of 
the Jovian atmosphere, and implies that observations in the 7.7| i CH^ band 
will provide relatively low vertical resolution information on the tempera- 
ture profile. 


2.2.2 10. 5u NH 3 Band 

For this band, we follow Gille and Lee (1969), who calculated 
the positions, half -widths and intensities of absorption lines due to 
ammonia at v <1400 cm \ and used this information to calculate the 
parameters of the Goody random absorption model with exponential distri- 
bution of intensities. For such a band, the transmittance can be written 
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where a is the absorber amount (cm STP) , P g is the effective pressure 
for line broadening (atm), and N is the number of lines in the spectral 
interval of width Av. It is assumed that the half width of the ith line 

i 

at temperature 6 is given by 


a.(e> = a 0 (e r ) f(e,e r ) P e 


( 22 ) 


where the subscript o refers to a pressure of one atmosphere of NH^, and 
f(0.0 r ) is a temperature correction factor. 


If we let 


N 

cr = y S.( Q ) 
i=l 


(23) 


and 


N 

£ '(•,.<» « 0i ) 
i=l x x ' 


1/2 


(24) 


we can write (19) as 


T = exp 


aa I" a 1 q 

Av [4 P e f(0,0 r ) p 


+ 1 


- 1/2 


(25) 


Based upon measurements of NH^ lines in the microwave region, 
Gille and Lee assume that 


P • t + P H / 7.9 + P / 11.7 
e NH 3 H 2 H e 


(26) 


For the temperature dependence, they adopt 
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(27) 


f(0 


’V ■ (r) 


5/6 


for hydrogen broadening. 


Based upon tjieir calculations, Gille and Lee tabulate the 
values of cr and p as a function of wavenumber from 12.50 cm ^ to 
1387.50 cm ^ in 25 cm ^ intervals, for temperatures between 100°K and 
225°K at 25°K intervals. We reproduce here as Table 2, the portion of 
their tabulation of application to the present study. 


We may use the strong line approximation of Equation (25) 


T = exp {-2 Av" 1 P £a P g f(0, 0 r )J 


1/2 


if 


D — Si 

~ 4 p f(e 
e 


I (zX 

(e.e r ) 


> > i 


For the Jovian atmosphere, we have 
a = 700 cm STP of NH 0 


P g = p atin /10 (for a predominantly E^, H g atmosphere) 

=0.2 atm near the cloud level 
5/6 . 


" 155 


= 2 


~ -1 2 
at v = 1000 cm , cr = 10 , p = 10 


For these values, we obtain 


R = 4.4 x 10 > > 1 


(28) 


(29) 
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TABLE 2 


VALUES OF cr AND 0 FOR THE 10. 5^ BAND OF NH 3 
(T = 150°K, Av = 25 cm *) (After Gille and Lee, 1969) 


V 

(cm V 

(microns) 

, -1 
(atm 

u 

cm -2 ) 

, „ -1/2 
(atm cm 

•3/2)* 

837.5 

11.9 

3.03 

X 

10 

1.33 

X 

10 

862.5 

11.6 

8.68 

X 

10 

1.81 

X 

10 

887.5 

11.3 

7.03 

X 

10 

1.12 

X 

10 

912.5 

11.0 

5.47 

X 

10 

1,20 

X 

io 

937.5 

10.7 

3.25 

X 

io 2 

5.91 

X 

10 

962.5 

10.4 

3.55 

X 

10 2 

6.62 

X 

10 

987.5 

10.1 

9.54 

X 

10 

9.74 



1012.5 

9.88 

1.37 

X 

io 2 

1.64 

X 

10 

1037.5 

9.64 

2.34 

X 

io 2 

3.19 

X 

10 

1062.5 

9.41 

9.93 

X 

10 

2.14 

X 

10 


There appears to be an error in Gille and Lee's (1969) Table 2b, 

-1 -3/2 

where the units of g are given as atm cm 
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Hence, we may use the strong line approximation of the 
Goody random model for the 10.5|j. NH^ band in applications to the Jovian 

atmosphere. With the substitutions Av = 25 cm \ P g = rP, where P = 

— 1/2 ~ 

P . , the average pressure along the absorbing path, and [f(6,0 )] = 

3Ctti r 

1.4, Equation (28) becomes 


[- 


T = exp -.112 £ (a r P) 


>1 


1/2 


(30) 


The amount of absorber and effective pressure of a vertical 


path of NHg can be written as 


00 

I 


“V 2 


a = 


STP 


(31) 


and 


00 

/ 


P da 


P = 


a(z) 


(32) 


Thus , 


CO 

- / 


aP = / P da 

z 


w 

-J- / 

STP J 


n NH_ 


dz 


(33) 


Assuming a scale height H for the distribution of pressure 
above the cloud-top (c), and a different scale height K for the distri- 
bution of NH^ above the cloud-top (since the NH^ is probably saturated 
in at least part of the atmosphere above the major cloud level, and, 
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hence, its number density is an exponential function of the temperature, 
with a different scale height than that of the atmosphere), we obtain 


aP 


p (c) (c) 

n STP 



(34) 


If Equation (34) is substituted into the formula for the 
transmittance, eq. (30), we obtain 


T = exp <( -A I J exp (-z/D) dz 
z 


1/2 


(35) 


where 


A = 0.112 p |r P(c) n^^ ( c > /n STp] 


1/2 


(36) 


and 


I-I + I 

D K H 


(37) 


Upon integration of (35), we obtain 


[• 


1/2 


T(z) = exp -AD exp (-z 


s/2D)] 


(38) 


from which the weighting function is derived as 


g - 0.5 Mf l/2 e' s/2D exp [-AD 1/2 


(39) 


For our preliminary calculation of the NH^ weighting functions 
in the Jovian atmosphere, we have assumed the following values for the 
Jovian atmospheric parameters: 
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n NH 3 (c) = 1.1 x 10 17 
a = 700 cm STP 

From above values , K = 1 . 72 km 
H = 1 6 . 5 km 
0/m = 51.5 
r = 0.1 

P(c) = 2 atm and 4 atm 

Weighting functions for three wavenumbers in the 10.5|i NH^ 
band are plotted in Figure 2. In comparison to the CH^ weighting functions, 
the NH^ weighting functions are much narrower and cover a smaller alti- 
tude range. This stems from the fact that due to condensation the NH^ 
scale height is much less than the CH^ scale height. 

It is of interest to see what the NH^ weighting functions 
would look like if the NH^ were uniformly mixed with altitude with the 
same scale height as CH^. In such a case, equation (33) becomes 

00 

a P = — — / P n dz (40) 

n STP * 
z 

where w is the NH^ mixing ratio equal to n^ /n, where n is the Jovian 
atmospheric number density. With the use of the hydrostatic equation, 
and the relationship 

dp = 10 6 dP (41) 
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Figure 2 . Weighting functions in the Jovian atmosphere for selected wave 
numbers of the 10.5(i NH^ absorption band. 



where p is atmospheric pressure in cgs units and P is atmospheric pres- 
sure in atmospheres, we obtain 


aP = 


io^w r 

n STP “ 


P(z) 

P dP 

o o 


(42) 


Integrating Equation (42), and substituting the result into 
Equation (30) for T, we obtain 


T = exp 


-0.112 x 10 3 p p' ^ 


feri 


(43) 


Letting the optical depth 


t = 0.112 x 10 J p. P 


( 2 "stp ) 


1/2 


(44) 


differentiating (43) with respect to height, and making use of the equa- 
tion of state, we obtain the weighting function 


dT 0.112 x 10 8 P ( rw tig 


) 


1/2 


- T 


(45) 


In evaluating (45) , we use a value of w compatible with the 
same total amount, (a), of NH^ es used in the previous calculation, that 
is, 700 cm STP. This value of w may be obtained from 


f 

STP W 


n NH dz 

*STP ^ wtl 3 
z 


CO 

_ w r 

- “ / n 

STP J 


dz 


(46) 
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With the use of the hydrostatic equation, (46) becomes 


From which, 





n STP 


( 47 ) 


w 


3 n STP 


(48) 


Taking p equal to 2 atm, we obtain w = 1.32 x 10 



The weighting function for the case of constant NH^ mixing 

ratio is compared to the weighting function for the case of decreasing 

NH^ mixing ratio in Figure 3. For the case of constant NH^ mixing ratio, 

the weighting function peaks at a much greater altitude (70 km versus 

-2 -1 

13 km), is much broader, and has a much lower peak value (2 x 10 km 
- 2-1 

versus 12 x 10 km ) . It is quite obvious that the region of the 
atmosphere to which 10.5pi emission observations refer depends critically 
on whether or not NH^ does indeed condense out and, hence, has a smaller 
scale height than the major atmospheric constituents. If, as is probably 
the case, NH^ does condense and follow a distribution with a small scale 
height, the narrow weighting functions that result will permit good 
vertical resolution in the deduction of temperatures from emission obser- 
vations. 


2.2.3 Pressure Induced Absorption 

For the pressure induced absorption of hydrogen, which, at 
Jovian atmospheric temperatures, peaks between 400 cm * and 800 cm ^ , 
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Figure 3 . A comparison of weighting functions for NH„ 

for two different assumptions on the distribu- 
tion of NH^ with altitude. 
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we follow the work of Trafton (1967). The volume absorption coefficient 


of at wavenumber v in a hydrogen -helium mixture can be written as 


A 


V , V c 


-1 


( A v + ’ B v) 


(49) 


where the units of ^ are cm , n^ is the H 2 number density, c is the 
speed of light, q is the number density ratio of helium to hydrogen, and 
A and are the pressure-induced binary absorption coefficients cor- 
responding to absorption in pure Hj and the enhancement in He-I^ mixtures, 
respectively. 


The transmittance may be written as 


T = exp 


’ CO 

A 


dz 


L z 


(50) 


where we have omitted the subscript v for convenience. 


Substituting the expression for k^ ^ into Equation (50), we 


obtain 


T = exp 


00 

- c / <A + -1 B) v 2 dz 


2 


(51) 


A and B vary womewhat with temperature, and hence, would vary 
with height. However, for an isothermal atmosphere, they are independent 
of height and Equation (51) becomes 
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T = 


exp 


. 1(A + q B ) 



dz 


(52) 


The number density of ^ above the clouds can be written as 


^2 = “h ^ 6XP 


(53) 


Integrating Equation (52) with the help of (53), we obtain 


T = exp 


(A + q B) ny (c) H exp(-2z/H) 
2c 


(54) 


From Equation (54), the weighting function can be derived as 


d| = lA_±_a_Bl 2 (c) eX p ( . 2z /H) e‘ T (55) 

where 

x = 2 (c) H exp(-2z/H) (56) 

We have computed the Hj weighting functions for the following 
values of the parameters. 

A and B obtained from graphical representation of 
Trafton (1967) for 9 = 160°K. 

q = 0.5 
H = 16.5 km 
6/m = 51.5 

P , , = 2 atm and 4 atm 

cloud 
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Shown in Figure 4 are the weighting functions for ^ for 
v = 700 cm ^ and v = 800 cm These wavenumbers are of interest in 
the analysis of the Jovian infrared spectral observations of Gillett 
et al. (1969), which cover this region of the absorption band. The 
large overlap of the two weighting functions is due to the fact that the 
absorption coefficient does not vary greatly in this spectral range. 

2.2.4 NH ^ Microwave Absorption 

For the microwave spectrum of NH^ there appears to be no 
simple analytical absorption model from which the weighting functions 
can be derived. However, information on absorption coefficients and 
half-widths for individual lines, and pressure broadening, is available 
in the literature (Townes and Schawlow, 1955; Anderson, 1950; Birnbaum 
et al., 1953; Legan et al., 1965; and Wrixon et al., 1971). To derive 
the transmittance at a particular wavenumber from this information re- 
quires an integration of the contributions from all the lines that in- 
fluence the absorption at that wavenumber. This is best done by an 
electronic computer and we plan to program such computations during the 
next quarterly period. 

2.2.5 Preliminary Estimate of Jovian Temperature Profile 

Assuming that the observed brightness temperature at a 
particular wavenumber represents the temperature at the altitude at which 
the weighting function for that wavenumber peaks, we can estimate the 
temperature profile in the Jovian atmosphere from the Jovian infrared 
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Figure 4- Weighting functions in the Jovian atmosphere for selected wavenumbers 
of the pressure induced H„ absorption band. 




spectrum (resolution AL/X = 0.02, or ~ 20 cm ^ at 8p.) observed by Gil- 
lett et al. (1969). We use the weighting functions computed above for 
CH^, NH^j and • Brightness temperatures were read from Hogan et al.'s 
(1969) plot of Gillett and Low's data. Table 3 lists the observed bright- 
ness temperatures and the altitudes to which they refer for selected 
wavenumbers of the three absorption bands. Two different height scales 
are listed: one relative to a cloud ai p = 2 atm and one relative to a 

cloud at p = 4 atm. For both height scales, a scale height of 16.5 km 
is assumed. 

Based upon these temperatures, a preliminary temperature pro- 
file for the Jovian atmosphere is plotted in Figure 5. The cloud tempera- 
ture of 210°K is based upon the results of Owen (1965; T c i ou( i = 200 + 25°K) ; 
Danielson (1966; = 210 + 15°K; and Gillett et al. (1969; T cloud = 

220°K) . A cloud top pressure of 4 atm is assumed in this plot. This 
yields a lapse-rate of in the first 20 km above the cloud ~ 4°C/km. 

The adiabatic lapse-rate for a 0.5 Hg and 0.5 He atmosphere is 3.3°C/km; 
for a 5/6 ^ and 1/6 He atmosphere, 2.4°C/km. Had we assumed a cloud 
pressure of 2 atm, we would have obtained a lapse rate of ~ 9°c/km in 
the first atm, far greater than the probable adiabatic lapse rate on 
Jupiter. Thus, the inferred temperatures appear to be consistent with 
a thick cloud at a pressure or 4 atm or more. 

The deduced Jovian temperature profile is similar in shape 
to the mean temperature profile in the earth's atmosphere. That is, a 
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TABLE 3 


OBSERVED JOVIAN BRIGHTNESS TEMPERATURES (AFTER GILLETT et al. 1969) 
AND THE ALTITUDES TO WHICH THEY REFER FOR SELECTED WAVE- 
NUMBERS IN THE CH, , NH , AND H„ INFRARED ABSORPTION 
BANDS J 


Gas 

V 

/ ~ l \ 
(cm ) 

t b 

(°K) 

^cloud 

z (km) 
4 atm 

for : 

P , = 2 atm 

cloud 

CH, 

4 

1180 

140 

41 


30 


1200 

138 

41 


30 


1225 

137 

51 


40 


1250 

140 

76 


65 


1300 

148 

91 


80 

nh 3 

875 

127 

20 


9 


1050 

125 

22 


11 

H 2 

800 

126 

22 


11 


700 

118 

26 


15 



HEIGHT ABOVE CLOUD (km) 


0I6A207-30E 



Figure 5 . Preliminary estimate of temperature profile 
in the Jovian atmosphere from analysis of 
infrared emission observations. 
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tropospheric layer, in which the temperature decreases with altitude and 
a stratospheric layer in which the temperature increases with altitude. 

The presence of such a Jovian stratospheric layer in which the tempera- 
ture increases with altitude was first suggested by Gillett et al. (1969), 
on the basis of a qualitative interpretation of their Jovian spectral 
observations. Theoretical calculations of the Jovian temperature. prof ile with 
a radiative-convective equilibrium model (Hogan et al., 1969) indicate 
that the observed stratospheric structure is due to absorption of solar 
infrared energy by the 3020 cm ^ band of methane, a suggestion made by 
Gillett et al. (1969). 
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3. FUTURE PLANS 

Plans for the next quarterly period include the following: 

(1) Programming of the computations of the NH^ microwave weighting 
functions for the Jovian atmosphere. 

(2) Further analysis of the Jovian temperature profile, utilizing 
available observations. 

(3) Commencement of development of inversion techniques for appli- 
cation to planetary atmospheres emission observations. 
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